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Abstract 

In this essay we examine the response of an Unruh-DeWitt detector (a quantum two-level system) 
to a gravitational wave background. The spectrum of the Unruh-Dewitt detector is of the same 
form as some scattering processes or three body decays such as muon-electron scattering or muon 
decay. Based on this similarity we propose that the Unruh-DeWitt detector response implies a 
“decay” or attenuation of gravitons, G, into photons, 7 , via 0 + 0 —>7 + 7 or G *—>7 + 7 +G*. 
Over large distances such a decay/attenuation may have consequences in regard to the detection 
of gravitational waves. 
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In the absence of an acceptable theory of quantum gravity the semi-classical theory of 


quantum fields in curved space-times 


-3] is the best available approach to the study of the 
interaction of quantum fields with gravitational phenomena. One of the key results of this 
semi-classical approach is the phenomenon of particle creation in a gravitational background. 
The usual example of this is Hawking radiation ( 4 ] - the creation of photons in the background 
of a black hole. One of the most direct ways to check if field quanta/photons are created by a 
given gravitational background is through the use of an Unruh-DeWitt (UD) detector 3, 5]. 
An UD detector is a quantum system with two energy levels. A simple, realistic model of 
an UD detector is an electron placed in a uniform magnetic fields [6J. By placing the UD 
detector in a given gravitational background, and observing how the upper energy level of 
the quantum system is populated, one can determine the emission/absorption spectrum for 
a given space-time. The UD detector method can be used to find the radiation associated 
with several well-known space-times: (i) the thermal radiation, with Hawking temperature 
T = 8 kg^MG ’ a black hole of mass M; (ii) the thermal radiation, with temperature 
T = 


ha 


27 TCki 


-, observed by an eternally accelerating Rindler observer. 


In this essay we want to examine the response of an UD detector to a gravitational plane 
wave background. The background of a gravitational plane wave is different from a black hole 
space-time in that the gravitational plane wave has no horizon. In this regard the response 
of an UD detector to a gravitational plane wave background is similar to the response of 
an UD detector undergoing circular motion as studied in (<3-9]. Neither the rotating space- 
time nor the gravitational wave space-time has a horizon (although the circular rotating 
metric has a light surface). The spectrum of the UD detector in the background of a plane 
gravitational wave is obtained and found to have the same form as some Standard Model 
scattering or three body decays. Based on this similarity we propose that this represents, in 
Feynman diagram language, either a “decay” of the initial two gravitons into two photons 
(be. G + G —* 7 + 7 ) or the attenuation of one initial graviton into two photons plus a 
lower frequency/energy graviton [i.e. G —* 7 + 7 + G). Such processes, as well as related 
processes like graviton-photon “Compton scattering” G + 7 —» G + 7 , have been investigated 


m 


10 


12] using using Feynman diagrammatic methods or other techniques. The attenuation 


process, G —» 7 + 7 + G, is related to the previously studied graviton decay [13j process 
G —» G + G + G. Both G + G —» 7 + 7 and G —* 7 + 7 + G would lead to a weakening of 
the gravitational wave amplitude/graviton intensity with distance from the source. 
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In order to examine particle production from quantum vacuum fluctuations we consider 
the special case of a gravitational plane wave [I14J traveling in the £3 = +z direction with 
diagonal metric. Taking the UD detector to be located as z = 0 the metric at this point is 
g 00 = —1, gu = 1 — h(ut), g 2 2 = 1 + h(ut), g 33 = 1 with h(ut) = sin 2 (0) sin (2'0) h 0 sin (ut). 
The angles 9 and ^ represent the orientation of the direction of the gravitational wave with 
respect to the detector axis. These angles take values from 9 = 0...7r and i/j = 0...2-7T jis| . 
As an example, for an UD detector consisting of a particle with a magnetic moment in a 
uniform magnetic field the angles of orientation, 0 and are between the direction of the 
magnetic held and the direction of travel of the incoming plane wave. The constant /i 0 is 
the dimensionless amplitude of the gravitational wave, which can be written as /i 0 = G 
where m 0 and r 0 are the relevant mass and distance of the system and is a potential 
energy per mass. The amplitude squared of the wave is proportional to the energy density 
jit] ]. Too — The graviton is taken as the quanta of the gravitational wave with 

energy hu. 

The spectrum, S(E), observed by the UD detector can be viewed as a production rate of 
photons per unit volume. It is given by the product of the density of states, p(E), and the 
response function F(E) 


(i) 


Tlgeneral inertial 27T /?(-£/) F (-£/) 

= 2 ME) /_+” e~ iAEA - (G + (s) - G+(Ar)) d( At). 

In this expression for S(E) we have taken h — 1 and c = 1. In ([T]) A E = E up — Edown is the 
difference between the energies of the upper and lower energy level of the UD detector jg]. 
We will take Edown = 0 so that from now on we can write A E —> E up = E. The number 
density of the inertial detector, rii ner uai , has been subtracted out to eliminate the singularity 


18] of equation (jT]) . The trajectory for a particle in an inertial frame 


in the integrand [7, 

is x(At) = (Ar,0,0,0) which has an associated Wightman function of G + (Ar) = yyryyr- 
The particle trajectory in the curved space of a passing gravitational wave is given by 
x(Ar) = (7 At, Ax, 0, 0) where y~ 2 = 1 — Ax 2 and Ax is the deviation of a particle’s 
trajectory due to the gravitational wave background. Substituting this into the general 
number spectra in equation (Op) the integral becomes 


3 


















( 2 ) 


F(E, 0 , iP) = /_ + “ e~ iEAr ( G + (s ) - G+(r)) d(Ar) 

= 2 ^ io°° cos (^ Ar ) ( 7 2 A ^-a,^ - a?) rf ( Ar )- 
Equation (J 2 j) is correct in general for strong gravitational waves since up to now we have 
placed no restrictions on ho or 7 . The dependence of F(E, 9, ip) on the angles 6 , ip comes 
from Ax as we will see shortly. Without any approximation this F(E,9,ip) for the plane 
wave background, as given in equation (J 2 ]) , must be evaluated numerically. 

In order to make an analytical evaluation of F(E, 9, ip) we take the limit of a weak 
gravitational wave with 7 « 1 and h 0 << 1. To first order the trajectory along a null 
geodesic [14] in the x direction is x = (l + \h). The separation, Ax , can be calculated from 
the difference in the null geodesic between particles in the inertial and curved space-times, 
Ax = (l + |/i) At — At = \hA r. Substituting this Ax into equation (J 2 ]) and using 7 « 1 
and ho << 1 we obtain 

F ( E > V0= 8^ r cos ( EAr ) Ar^t- 1 / l 2 ) d ( AT ) ~ 8^ Jo°° COS( fy fe2 d(Ar). (3) 

With this approximation the RHS of equation (J3]) can be integrated. Recalling that 
h(AT,9,ip) = sin 2 (9) sin (2ip) ho sin (u At) (this dependence of h on 9, ip is the reason F 
in (J2J and (J3J) depends on 9, ip) we find, 

Finally we take the integral of (J3J) over all possible orientations of the direction of the 
gravitational wave with the orientation of the UD detector, (he. integrating over 9 = 0...7T 
and ip = 0...27r) to find 


y sin 4 (9) sin 2 (2ip) h^. (4) 


F(E) = ^ (2uj — E) hi for E < 2hu ; and F(E) = 0 for E > 2hu>. (5) 


The response function F(E) vanishes for energies above the cutoff E > 2huo. Note that this 
cut-off is also the same as that seen in muon decay 19] 20]. The meaning of the cut-off in 
(|5]1 can be seen to have a similar kinematic origin - for the processes G + G —> 7 + 7 and 
G —» G + 7 + 7 there is at most 2 hu of energy available (each graviton contributes hu of 
energy). Thus for E > 2 hu neither process can proceed energetically. 
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The density of states in equation (JT]) is p(E) — ^ [7|. Using this along with the response 
function in equation (J5J), and restoring the factors of h and c, gives the spectra detected by 
the UD detector 


S{E) = 


5127r/i 3 c 3 


E 2 (2kuj — E) hi for E<2hu ; and S(E) — 0 for E>2hu. ( 6 ) 


The result in (J6J) gives the spectra of the UD detector S(E). Equation (J6]) is the main result 
we want to emphasize in this essay. As previously mentioned S(E) has the same form as 
certain Standard Model processes (see for example 19] 20]). The functional form of S(E) 
in equation ([ 6 ]) is a Beta( 3,2) distribution and such distributions are common for particle 
decay processes. In Feynman diagram language, the processes implied by S(E ) from (J6]) 
would beG + G —)-7 + 7 orG—)-G + 7 + 7 where the gravitons take the place of the muon 
and electron and the two photons take the places of the neutrinos. For G + G —)■ 7 + 7 
the gravitons are completely transformed into photons while for G —> G + 7 + 7 the initial 
graviton is transformed into a lower energy graviton. For either process energy is transferred 
completely or partly into photons and the gravitational wave will be attenuated. 

Recalling that S(E) from (J 6 ]) can be viewed as a production rate of photons per unit 
volume for the UD detector we can use it to calculate the decay/attenuation of gravitational 
waves via the above processes. By energy conservation the photon production spectrum 
for the UD detector will come from the gravitons. The “decay” rate for gravitons per unit 
volume is then given by, 


2 huj 

= J S(E)dE. (7) 

0 

The production rate is a function of the dimensionless gravitational wave amplitude, 
h 0l but the physical amplitude falls off as r -1 for an expanding gravitational wave. The 
dimensionless gravitational wave amplitude can be written as h 0 = G where mo and r 0 
are an effective mass and distance respectively. The physical amplitude for an expanding 
gravitational wave can then be written as h r = y where £0 = TprU which has dimensions of 
length. Replacing h 0 with h r in equation ((7J) and integrating over all energies in the spectra 
yields, 
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dT = uk?> idV = < 8 > 

One can immediately integrate the RHS of equation (JHJ) from R = 0 to some radius R = r 

c2 

to obtain a decay rate of T = 3 ^ 3 -co r. The four factors of h coming from hw are canceled 
by the three factors of h coming from S(E) and the one factor coming from the LHS of 
equation ((7|). 

The decay of gravitons can be expressed in terms of a mean path length, (r) = c (r) = 
cT” 1 . Integrating equation (JHJ) from 0 to (r) and substituting T = in the mean path 

length expression gives, 


. . c 32c 4 
<r> “ f “ (f? (r) 


( r ) = 


2a/8c 2 


( 9 ) 


Due to the inverse relation between the mean free path, (r), and £ 0 an d w 2 , gravitational 

waves are effected less, and therefore attenuated less, by this conversion of gravitons into 

photons for low frequency (u <C 1) and small amplitudes (£0 1). This is similar to 

the conclusion reached by Mondanese [ 13 I using “general kinematical arguments”. For a 

gravitational wave source on the order of a Solar mass, mo ~ 2 x 10 30 kg, we find £0 ~ 10 3 m. 

Next assuming a gravitational wave frequency of co ~ 10 2 s _1 , the approximate mean decay 

length would be (r) ~ 10 n m = 10 ~ 5 ly. This is a short attenuation distance especially 

considering that the decay rate, T, is normally relatively low. However, this short attenuation 

length is not completely surprising since a small decay rate can be compensated for if it acts 

over a large distance traveled by the gravitational wave. Also for strong gravitational wave 

£ 2 

sources with large amplitudes, £ 0 > an d large frequencies, uj, the decay rate, T = 
becomes larger. Thus for strong gravitational waves one would need to numerically integrate 
()2|) without the approximations and then repeat the above development. 

The mean free path for graviton decay in equation (jUJ) is based on a semi-classical calcula¬ 
tion with the gravitational wave treated classically and the vacuum quantum mechanically. 
The quantum nature for the graviton appears here only via the particle energy, hui, and 
through the cutoff in the energy distribution, 2hu. It would be interesting to compare 


this resu 
DeWitt 


t to a more complete quantum held theory calculation of graviton decay following 


2l| or Bjerrum-Bohr et al. 12]. While such calculations have the advantage of 


including spin and more carefully taking account of energy-momentum conservation, any 
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mean free path calculation for graviton transition rates would include a similar calculation 
over the volume of space as in equation (JSJ). The UD detector calculation, followed by the 
mean free path estimate given above, indicates that detecting gravitational waves at large 
distances may be hard or impossible due to the attenuation of the gravitational wave coming 
from the conversion of gravitons to photons. 
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